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Abstract 

In this paper, we derive a simple sum rule satisfied by the gluon spec- 
tral function at finite temperature. This sum rule is useful in order to 
calculate exactly some integrals that appear frequently in the photon or 
dilepton production rate by a quark gluon plasma. Using this sum rule, 
we rederive simply some known results and obtain some new results that 
would be extremely difficult to justify otherwise. In particular, we derive 
an exact expression for the collision integral that appears in the calcula- 
tion of the Landau-Pomeranchuk-Migdal effect. 

LAPTH-909/02, LPT-ORSAY-02/27 



1 Introduction 

Photon production is considered to be a very interesting signal of the formation 
of a quark gluon plasma in heavy ion collisions Indeed, because 

of their very weak coupling to matter, photons (and more generally any elec- 
tromagnetic probe) have a large mean free path which enables them to escape 
without reinteractions from a medium the size of which is at best a few tens of 
fermis. 

On the theoretical side, the calculation of the photon and dilepton rate is 
performed under the hypothesis of local equilibrium, i.e. one calculates a rate 
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Figure 1: The two- loop diagrams contributing to photon and dilepton produc- 
tion. 





Figure 2: Important processes contained in the above 2-loop diagrams. 

(number of photons produced per unit time and per unit volume) using thermal 
field theory in equilibrium, and then plugs this rate into some hydrodynamical 
model which describes the system by dividing it in cells where a local 

equilibrium is assumed and by assigning a local temperature and fluid 4- velocity 
to each cell. Such a description is consistent as long as the photon formation 
time is small compared to the typical size of the cells in which an approximate 
equihbrium is reahzed 

Thermal field theory calculations of photon and dilepton production rates 
have been performed a long time ago [[lOyil| , |l^ , p^ , and have been reassessed 
under the new light shaded by the concept of hard thermal loops (HTL - 
|l5|,^6ljl|,|l|l) |l|,||,|l|j2|,|2|,|j. In this context, it has been found that some 
2^3 and 3^2 processes which appear only in 2-loop diagrams [p5|j2^ , |27| , p8| 
(like bremsstrahlung, as well as the process qq{q^ g} — > 7{'7, g} which can be 
deduced from bremsstrahlung by crossing symmetry - see figures and |) are 
also important. These processes are in fact enhanced by a strong sensitivity 
to the forward emission of the photon, due to a collinear singularity regularized 
by a thermal mass of order gT where g is the gauge coupling. In particular, 
the process on the left of figure || has been found to enhance considerably the 
rate of hard photons pTf . This collinear enhancement mechanism has also been 
shown to play a role in multi-loop diagrams belonging to the class of ladder 
corrections and self-energy-corrections [ p9| , ^0pl| ,|8|, p2||33| , and a resummation 
of this family of diagrams has been carried out in p2| , |33| . The effect of this 
resummation, also known as Landau-Pomeranchuk-Migdal (LPM) effect p4,p5l 
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p6[ , leads to a small reduction (by about 25% for photons in the range interesting 
for phenomenology) of the photon rate. 

However, the strength of the LPM suppression decreases as one increases 
the invariant mass of the photon since the photon mass helps to regularize the 
coUinear singularities. It is therefore expected that the dilepton rate can be ac- 
counted for by limiting oneself to a 2-loop calculation, and that the production 
of hard dileptons of intermediate invariant mass is dominated by the process 
shown on the left of figure ^. Under this assumption, the calculation of |^ 
has been extended recently to the case where the photon mass cannot be ne- 
glected any longer, with emphasis on the "off-shell annihilation" process that 
was already dominant for hard photons ||3^ . In this paper, we found a sim- 
ple generalization to the case of massive photons of the formula known for the 
imaginary part of the real photon polarization tensor. This formula reads: 

Imn^^(g)w • — / dpo[n^iro) -n^ipo)] 



2n^ go' 



— OO 



with ro = po + qo and^ 



1 -l-oo 



dx^ ^ , , f , \/wl{w + 4)tanh ^Jwl(w + 4) 

I t-v-\ II / m 1 f y-i rt i i * * 



/dx f 

— Im , (x) / dw , , , , > , o . 

X ' J (M^V + Ren, Ja;))2 + (Imn,,(x))2 ' 

' 

1 +00 _ 

-Af2 fdx. „ , . [dw ^w/{w + A)t&nh ^ / {w + A) - w / A 







w (M4u;-KRen^^(x))2 + (Imn^^(a;))2 ' 

(2) 



where the 11^ ^ are the transverse and longitudinal self-energies resummed on 
the gluon propagator, and where we denote 

M^fj EE Ml + ^poro , (3) 
% 

with Moo the thermal mass of a hard quark {M^ — g^CpT^ /A). The functions 
^ are the usual transverse and longitudinal HTL gluon self-energies: 



'"x^ x{l 



In 



x-l 



^Let us recall here that M^jj can become negative if > 4M^. Those definitions are 
only appropriate for M^g > 0, because they have been obtained by rescaling thji_transverse 
momentum of the gluon by writing fi^ = M'^^w. However, it was argued in |37| that the 
correct result for M'^^ < can be obtained as the real part of the analytic continuation of the 
result for > 0. Most of this paper deals with the case of a positive M^^. 
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(L) = 3ml 

(4) 

where we denote x = Iq/I and where mg = g^T^[Nc + Np/2]/9 is the gluon 
thermal mass in a SU{Nc) gauge theory with flavors. 

Therefore, in addition to the function J^^ already introduced in the case of 
quasi-real photon^ we need two new functions K^^ for the term proportional 
to the photon invariant mass squared Q^. All are dimensionless functions of 
the ratio of M^s to the plasmon mass rUg which appears as a prefactor in the 
self-energies ^ . Up to now, the ^ and , have only been evaluated 
numerically ( ||26| , p7| , with a mistake corrected by [ p8l|39| , ^ ) , which is sufficient 
for the case of real photons since in this case they are fixed numbers that depend 
only on the number of colors and flavors but not on kinematical parameters 
(M^jj = M^). However, the cost of this procedure increases significantly in the 
case of virtual photons since the value of M^g depends on the invariant mass 
Q, energy go a,nd quark energy po- In addition, obtaining asymptotic limits is 
far from trivial at this point]^. 

The aim of the present paper is to derive some analytical results regarding 
those functions. We first show how the integral over the variable x can be 
performed exactly in the functions ^ and ^ by means of a simple sum-rule 
(section]^). This leads to either a very simple integral representation of these 
functions or even to closed formulas in terms of dilogarithms (section^). Thanks 
to these results, we can easily study the asymptotic properties of the functions 

^ and ^ . Non trivial asymptotic expansions are obtained, which would 
have been very difficult to obtain otherwise (section ^). In section || we show 
that the above analytic results can also give some insight on the fact that the 
processes of figure ^ depend only on parameters like the gluon screening masses 
and the hard quark thermal mass, in a generic model where the quark gluon 
plasma is described as a gas of quasi-particles. Finally, we show that one can 
also calculate analytically the collision integral that appears in the resummation 
of ladder diagrams |^,|3| (section^). In appendix ^ we derive the asymptotic 
behavior of a function introduced at an intermediate stage. In appendix we 
prove analytically an anecdotical property which was first noticed numerically 
the integrals and are exactly opposite if Moo — Wg and — 0. 



■^The term in M^K^, ^ was forgotten in M]. It eo mes fro m the HTL correction to the 
755 vertex. This vertex correction was also neglected in |32| , |33[ , without any damage to this 
approach since it affects only the component II^^ of the polarization tensor, while only the 
transverse components are calculated in these papers (see foij_more details on this issue). 

''Some very partial asymptotic results have been obtained in [E6f for j^. 
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2 Derivation of the sum-rule 



We want to calculate the integral: 

^ f' dx 2Imn(x) 
~ Jo X {z + Ren(x))2 + {lmU{x)y ' ^ ' 

for a positive z, where n(a;) is some self-energy depending only on a; = ko/k as 
is the case for instance with the HTL gluonic self-energy. The factor 1 /x comes 
from a Bose-Einstein factor in the soft approximation dkQn^{kQ) « Tdko/ko = 
Tdx/x. The first step in this calculation is to rewrite it as 

, dx ,^ 2Imn(a;) 

f(z)= — (l-a;^) _ ^ ' = , (6) 

' Jo x"- ^(z(a;2-l)-Ren(a;))2 + (Imn(x))2 ' 

where we define n(a:;) = (1 — x'^)Il{x). Interpreting now z (z > 0) as the square 
of a three-momentum k and x as the ratio x = ko/k, we have 

-2Imn(.x) 



(z(a;2 _ 1) -Ren(a;))2 + (Inin(a;))2 ~ k^ ~ U{ko,k)+ikoe 

= nko,k). (7) 

Note that this function is nothing but the spectral function p(fcoi k) associated 
with the "propagator" appearing in the right hand side. This is what enables 
us to relate the integral 

1 

f{z) = -J^{l-x')piV^x,V^) (8) 



to the spectral representation of this propagator. Indeed, it is known that the 
resummed propagator i/ {K^ — U{K)) is related to its spectral function p(fco, k) 
via the following spectral representation p2|j4^ : 

- / Ep(E,k)--^ 1 — . (9) 

k^-k^- U{ko/k) + ikoe Jo '^^ ' 'kl-E^ + ikoe ^ ' 

Taking the real part of this identity^, one recovers the definition Eq. (Q) of the 
spectral function. Taking its imaginary part and denoting E = kx and fep = ky, 
one obtains the following non-trivial integral: 



dE 



^xp{x, )y2^^2 (fc2(y2_i)_Ren(j/))2 + (Imn(?/))2 
Taking the limit y oo, we find 



do; _ 1 1 

—xpjkx, fc) - ^2 _ lim^^^ RcTi{y)/y^ ~ k^ + Ren(oo) ' ^^^^ 



Note that the spectral function p(ka,k) is by definition a real function. 
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Having in mind that !!(?/) is a gluonic self-energy obtained from the hard thermal 
loop approximation, its imaginary part vanishes if y > 1 and therefore does not 
contribute in the limit y oo. Alternatively, taking the limit y — > and 
assuming similarly that lmll(jj = 0) = 0, we obtain: 

'^°°dxl 1 

p(kx, k) = — . (12) 

TT x^^ ' ' P +Ren(0) ^ ' 

We can combine these two relations into 

'■+°° dx 







X 



{1 — X )p{kx, k) — TT 



fc2+Ren(0) /c2+Rcn(cx)) 



(13) 



In order to obtain from there the function /(z), we need to subtract the 
contribution coming from x between 1 and +oo. Fortunately, since Imn(a;) = 
for X > I, the contribution in this range comes only from the poles of the 
propagator i/iK^ - n(X)), via the formulaP 

^'-^il-Xrpikx,k)=TT Y: (14) 

poles Xi * 

where Z{xi) is the residue of the propagator at the corresponding pole. For the 
above propagator, the equation that determines the poles is 

fc2(a;2 - 1) = Ren(x) = (1 - x^)ReIl{x) . (15) 

This "dispersion equation" has a trivial solution x ^ 1, which does not 
contribute when plugged in Eq. (|lj) because the other factors in the integrand 
vanish if a; = 1. Any non trivial pole would be a solution of the equation 

k^ + Rell{x) ^0 , (16) 

but under the reasonable assumption that the resummation of the self-energy 
Il{x) leads to well behaved quasi-particles (i.e. that the equation kg — k^ = 
n(fco / k) has a solution for every value ofko/k larger than 1) , we have Re Il{x) > 
for a; > 1 and therefore there are no additional poles. As a consequence, the 
integral of Eq. ( |l^ ) does not receive any contribution from the range x G [1, +oo], 
and we can write directly a closed expression for the hmction f{z): 



^dx 21m n(a;) 



1 1 



z + Ren(oo) z-l-Ren(O) 



/o X {z + ReU{x))^ + (lmU{xyf 

(17) 

This is the basic sum-rule from which we are going to derive some results re- 
garding photon production by a quark-gluon plasma. The validity of this result 
can also be checked numerically. 

It may be useful to recall that, even if the derivation has been made having 
in mind a hard thermal loop for the self-energy 11, this result has a broader 
range of validity. In fact, it is valid for any self-energy satisfying the following 
assumptions: 

^If Imn(x) = 0, then p(fca;, k) = 2TTS{k'^{x^ - 1) - n(a;)). 
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1. n depends only on a; = k^/k 

2. Imn(2: = 0) = 

3. Imn(a;) = if x > 1 

4. Ren(x) > if x > 1. 

Note that the condition (1) can in fact be relaxed since what is done here 
can be reproduced if the self-energy depends separately on k and x, the only 
difference being that the result would depend on k. The condition (2) is always 
true. Conditions (3) and (4) depend on the nature of the resummation under 
consideration, but are reasonable approximations in any system of well defined 
quasiparticles. 

In p6) , a formula for the photon rate based on sum rules was also presented. 
However, in this paper, the use of sum rules led only to a very complicated 
result (involving explicitly the gluon dispersion equations as well as the residues 
of the gluon poles) that was not useful for any practical purpose. By comparing 
the two methods, we can trace the simplification achieved in the present paper 
to a different choice for the integration variables. Indeed, in ||2^ the sum rules 
were applied to an integral over the variables x = Iq/I, I = |'| (where L is the 
momentum of the exchanged gluon), while in the present approach, we take 
as independent integration variables x and w = —L^/M^g ~ P{1 — x^)/M^g. 
It appears that trading I in favor of w before using sum rules to perform the 
integral over a; leads to a dramatic simplification of the result, because some non 
trivial parts of the x dependence get absorbed in the new variable w. There are 
in fact many sum rules satisfied by the HTL spectral functions. The interested 
reader may find other examples in Q,^,^ . 



3 Expression of J^^ and 



T,L 



Thanks to the formula derived in the previous section, one can first simplify the 
functions and Kj, ^ so that we have only one-dimensional integrals over the 
variable w: 



+ C30 

— / dw 



w 



w -I- 4 



-tanh ^ 



w 



w + A 



1 



Ren (oo) 



A/2 



Ren^ ^(0) 



+ 00 



K^r =- — 



dw 
w 



w -I- 4 
1 



-tanh 



w w 



w + A 4 
1 

, Ron^^(o) 



(18) 
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At this point, using the change of variable u = yjw/{w + 4), we can write 
in a simpler way the following elementary integrals: 



+ OC 



W 







dw 

dw 

w 



where we define the function 

F{x 



-tanh 



w 



1 



w + a 



2F 



-tanh 



w 
1 



■In 



du- 



-F 



tanh ^ (u) 



(19) 



(20) 



{x - + 1 ■ 



Recalling now the following properties of the HTL self-energy of the gluon 

Re Uj. ^ (oo) = TOg 
Ren^(O) = 

Ren^(0) = 3TO^ (21) 

where mg is the plasmon mass, we easily obtain the following expressions in 
terms of the function F{x): 



= -ttF 

= TT 
Kj. — TT 
K=TT 



V mg / 



( m?- ) 



ml \ ml / 4 8 V ml ) 



3ln(3) 



F 



-If 

3 



3TOg 



(22) 



Therefore, those results demonstrate that in order to study the properties of 
the functions ^ and ^ , one needs only to study the properties of the much 
simpler function F{x). In fact, it is even possible to write the function F{x) 
in closed form in terms of dilogarithms^. We are not going to make use of this 



^Explicitly, we have: 
1 
4ip 



Li2 



+21n{2)ln 



2p 
p + i 

i — p 
i + p 



2Li2 



P 

p + i 



Li2 



2p 

p — i 



-2Li2 



with p= Vx — 1 and Li2(a;') = 



(23) 
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possibility here since it is simpler to keep F(x) in its integral form given by 
Eq. (|0|). 

One must stress the fact that all these functions depend only on the ratio 
of two masses. In the case of real photons {Q^ — 0), we have in addition 



M^g = and for Nc = 3 colors, we can write: 

4M2 8 



3ml 6 + 



(24) 



i.e. the temperature and strong coupling constant also drop out of this ratio. 
It appears that there is an additional (and purely accidental) simplification for 
Np = 2 flavors: in this case, the differences — Jj. and — can be 
expressed in a very simple fashion: 



K, -K. 



= 7rln(2) 



-(l-21n(2)) 



(25) 



For the case of = 3 flavors, the results are less simple, but one can still 
obtain explicit expressions: 

-2 33 



---ln^(2) 
8 8 ^ ^ 



3 ln(2) ln(3) 



3^. /3^ 
-Li 



2 

1 _ 7r2 
4 36 



:Li2 



ln(2) ln(3) 11 



8 



--In(2)ln(3) + -Li2 - 



— In (2) 
12 ^ ' 



1 

— Lio 
3 



(26) 



Had these exact formulas been known, the confusion due to the erroneous factor 
4 in the numerical evaluation of these coefficients in |^ would have been avoided 

m. 



4 Asymptotic behavior of J^^ and K^^ 
4.1 Limit Mcs ^ rUg 

Using the above results, we can recover in a rather simple and elegant way all 
the asymptotic limits given in ||2^ for Jj. j^, as well as the limits used for ^ 
in order to obtain the behavior of 2-loop dilepton production near the threshold 

= 4M^ (a region which is dominated by small values of M^g) \ ^7\ . 

To that effect, we need only the behavior of F{x) when x ^ 0. This is 
derived in the appendix where we prove that: 
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Thanks to this formula, a trivial calculation give; 




These relations in fact go well beyond the results for ^ obtained in |g6[ , since 
in this new approach we obtain for free the prefactor of the leading term and we 
could even have calculated some subleading terms, down to the constant term. 



4.2 Limit M^^ > 

Using Eqs. (p2|), it is also very simple to obtain the behavior of the functions 
^ and ^ in the opposite limit where A/cH 3> mg. In order to do that, we 
need to know the behavior of F(x) for large values of x. The following formula 
is also derived in appendix |^: 

From this formula, it is easy to obtain 

TT TO? . 



off 



g 



TT ml , /Ml 



J « -—^Infl^'l w -2J 
^ Mrft>mg 4 M^jj \ ml ) ^ 

K « ^^inf^ 

^ Meff>mg 48 A/2jj V m2 

TT ml / A/^£t \ 
« -;T7T7f In(^) « -2if^ . 30) 
Meff»mg 24 Af2^ V m2 / ^ ^ 



5 Beyond the HTL approximation 

In section |[ we mentioned the fact that the sum rule in Eq. is in fact valid 
for a gluon self-energy 11 more general than the standard case of hard thermal 
loops. Assuming it can be applied, we see that the result depends only on the 

^These leading log formulas are not affected by the fact that the asymptotic expansion of 
F{x) is slightly modified if x approaches with negative values (see Eq. ([79)). 
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four numbers Re 11^^(00) and Re 11^^(0). The first two are the plasmon mass 
(longitudinal) and the mass of the transverse gluon at zero momentum, and can 
be shown to be equal thanks to Slavnor- Taylor identities |^,^ 49 1 . Physically, 
this property means that there is no way to distinguish transverse and longi- 
tudinal modes for a particle at rest. Therefore, we need only to introduce one 
plasmon mass: 



Ren^(oo) = Ren^(cx)) 



(31) 



The quantities on the second line are squares of the screening masses for the 
transverse and longitudinal static gluon exchanges. The longitudinal screening 
mass is the familiar Debye mass: 



Ren JO) 



(32) 



In the HTL approximation, there is no screening for the transverse static gluons, 
but this is not expected to hold generally. The corresponding screening mass is 
the magnetic mass, and is denoted 



Ren^(O) 



(33) 



In terms of those parameters, it is straightforward to write down the ex- 
pressions of ^ and ^ for photon production in a description where we use 
gluon propagators that are more general than the HTL propagators: 



— -K 



J, = TT 



mag 
2 



V ml 



F 



V ml ) 



Af2 



F 



/r2 
off 



\ m. 



M3 



fF( 



AM. 



off 



\ m. 



(34) 



It is easy to check that these relations fall back to Eqs. (|2^) if we set mmag = 0, 
mp — nig and rriD — VSnig , which are the relations between masses in the HTL 
framework. 

There is another general property of the processes of figure |^ which is worth 
mentioning here. Their rate in fact depends only on the combinations J,, — Jj^ 
andif^-if^ (seeEq. @)) after one has summed the contributions of transverse 
and longitudinal gluons. Using the above formulas, we obtain: 



f( 



V mi 



8 \m^ J 



cff 



F 



"mag 



\m 



mag 



(35) 
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In other words, all the dependence on the plasmon mass drops out for the 
processes of figure ||. This property is in fact reasonable since we are looking at 
processes that involve only space-hke gluons, and it would have been surprising 
if the result had depended on the plasmon mass, a property of time-like gluons. 
The practical consequence of this for a phenomenological approach to photon 
production based on some quasi-particle picture is that we do not need to know 
the full gluon propagator, but only the two screening masses and the quark 
thermal mas^ Note also that these quantities remain finite even if the magnetic 
mass is very small or vanishing. 



In particular, it is now known from lattice calculations [pO 51 that the masses 



of quasi-particles increase when the temperature approaches the critical tem- 
perature from above, while at the same time the screening masses decrease. 
The formulas of this section are important to deal with such a situation, since 
they do not assume any particular relationship between the screening masses 
and the quasi-particle masses. For instance, for a temperature just above Tc, 
we can make use of Eqs. (^0|), and predict a suppression of the production rate 
of photons and low-mass dileptons simply due to the fact that screening masses 
are much smaller than the quasi-particle masses (in addition to the standard 
suppression due to the fact that the temperature is smaller). 



6 Resummation of ladder diagrams and LPM 
effect 



6.1 Integral equation in momentum space 



The authors of [ |32y33| ] performed the resummation of all the ladder diagrams, 
as well as all the self-energy corrections that are required to preserve gauge 
invariance, in order to account for the LPM effect in the production of real 
photons. Indeed, such photons may have a formation time larger than the mean 
free path of the quarks in the medium ]3l| |, so that multiple quark scatterings 
contribute coherently to the formation of the photon. 

In the formulation of [ p2[j3l^ , the imaginary part of the retarded photon 
polarization tensor is given byrT 



Imn^CQ) 



Stt 
xRe 



-f-OC 



Pi 



Pq^o 



(2^)2 



P± ■ f{Pi 



(36) 

*This is not true for the 2 — > 2 processes calculated in Indeed, since these processes 

involve time-like gluons, they can depend on the plasmon mass. 

^Note that this equation includes only the two transverse modes of the photon, and is 
therefore incomplete for the production of virtual photons. It is however not difficult to 
include the longitudinal mode as well l52|. 
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where = po + qa and f{p±) is a dimensionless transverse vector that repre- 
sents the resumnied couphng of a quark to the transverse modes of the photon, 
satisfying the foUowing integral equation 

2p^ = zSEf{p^)+g^C,T J ^C{l^)[f{p^) - fip^+l^)] , (37) 

where SE = qo{p± + Al^g) / {2poro) and in which C{1±) is the colhsion integral 
defined by 



2Inina(L) 



E 



=L,T 



(L2_Rcn«(i))2 + (inin„(L))2 



pnm^^Q^, (38) 



with Qf^ = {l,q/q) and P^'^^{L) the transverse or longitudinal projector for a 
gluon of momentum L. Note that SE~^ is nothing but the typical formation 
time of the photon |Q . Note also that in an iterative solution of these integral 
equations, the first term that contributes to the imaginary part of the photon 
polarization tensor is of order since f{p±) and g{p±) are purely imaginary 
at the order g^. This reflects the fact that the direct production of a photon by 
the processes —> 7 or g — + 57 is kinematically forbidden. This remark ceases 
to be valid if SE can vanish, in which case an ie prescription must be used, so 
that the zeroth order solution can have a real part. This happens if M^g can 
become negative, i.e. if > 4M^. 
Using the fact that we havepl 



(41) 



and introducing the variable x = Iq/I, we can rewrite^ the collision integral in 
Eq. (|8l) as 



2 r dx 

Jo X 



Im 11^^(2:) 



(/i+Renjx))2 + (Imn,(x))2 

Im Ily (x) 
"(^2^ + Ren^(x))2 + (Imn^(a;))2 



One can use 



and 



Q^Q,^ g' 



L2 



L2 ' 
if lo= Iz 



in order to obtain the second contraction. 

^^This change of coordinates has the following Jacobian: 



1 - 



'0 



din , o , dx , n , 

d(ll) = —dill) . 



(43) 

(39) 
(40) 

(42) 
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At this point, the sum rule derived in section g gives directly the result of the 
integral over the variable x, so that the collision integral can be rewritten as: 

We observe again that summing over the contributions of transverse and longi- 
tudinal giuon exchanges cancels the terms involving the plasmon mass. 

6.2 Solution in the Bethe-Heitler regime 

As a check, one can solve this integral equation iteratively up to the order 
g^. Indeed, the naive term of order is purely imaginary, and drops out of 
the imaginary part of the photon polarization tensor (see Eq. (|3^)). For the 
function f{pi_), this expansion gives: 

T3 f '^^P^ fl \ o -in rj.iPoro)'^ f (PPI- P 



Pi + Mis iP±+l±?+Mlff 



(45) 



At this stage, performing the angular integrations is elementary. Making use of 
the following identity: 

this can be rewritten as: 

d^P± . 2g^C,T {poroy 

(27r)2 0(g2) T,3 q2 

Finally, plugging Eqs. (^) into Eq. (|3^) gives 



-i-oo 

Imn,{:(g)^=^ ^^^7 y dpoK(ro)-n,(po)] 



— oo 



x{pl+rl)[J,-J,+2K,-2K,], (48) 

which is equivalent to Eq. @j for real photons {Q^ = 0, M^g = M^). This 
proves the agreement between the perturbative approach followed in p7[ and 
the resummation of the LPM corrections, if one formally keeps only the 0{g^) 
terms. The fact that some terms proportional to in Eq. ( p]) a rc not recovered 
in this limit is due to the fact that the LPM resummation of [ |32p^ is limited to 
the transverse modes of the produced photon, while massive photons also have 
a physical longitudinal mode. 
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6.3 Reformulation as a differential equation 

Since the collision integral appearing under the integral over l± is now known 
in closed form, it is possible to transform this integral equation into an ordinary 
differential equation by going to impact parameter space. We can first define 



Note that the order zero (in g^) solution of the integral equation is 
n qo TT go 



(49) 



(50) 



where the Ki's are modified Bessel functions of the second kind. For the higher 
order terms gi{b), one obtains the following equation for gi 



i^{Mls - AMb)+9'C,TD{msb){go{b)+g,{b)) = 



with 



D{m,b) 



2tt 



7 + In 



\/3ni„b 



ifo(\/3mg6) 



(51) 



(52) 



In addition, we have also 



Re p± 



cPpi 



■fiPA 



lim Im • gi(b) . 

6^0+ 



(53) 



This differential equation can be further simplified by defining the following 
dimensionless quantities: 



t ^ Ml^b^ , 

u{t)- 



b g{b) . 



(54) 



This transformation leads to 



Atu'{{t) - u^{t) + ig^C,?^^D i^Vt] {uo{t) + u,{t)) = , (55) 

where the prime denotes the differentiation with respect to t, and where 

uo{t) = VtKiiVt) . (56) 

The differential equation Eq. (^5|) depends on two dimensionless quantities. One 
is the ratio mg/Mcs of two masses, while the prefactor g'^C^.Tporo/qoM^g can 
be interpreted (up to logarithms) as the ratio of the photon formation time to 
the quark mean free path. It is therefore the average number of scatterings that 
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can contribute coherently to the production of a photon. The relevant quantity 
for the photon production rate is then given by 

This differential equation must be supplemented by boundary conditions in 
order to define uniquely the solution. First, we have 

u(0) = , (58) 

which can be seen from its definition. Subtracting uo(0) = Ij this implies 

Mi(0) = -1. (59) 

Unfortunately, we do not know the value of u'i{0) but instead we know the value 
of u{oo). Indeed, assuming that f{p±) is regular enough, its Fourier transform 
is exponentially suppressed at large b. Therefore, u(oo) — 0. Since we already 
have uo(oo) = 0, this implies 

ui{oo) = . (60) 

Therefore, the problem can be summarized as follows: we are looking for the 
(presumably complex) value of Ui(0) that takes us from ui(0) = — 1 to ui(oo) = 
0. The imaginary part of this derivative term is then the coefficient that enters 
in the photon polarization tensor. Note also that since the term and ui come 
with opposite signs in the combination Atu'( — ui in Eq. (p5|), generic solutionsf^ 
are unstable and diverge as i ^ +oo (they are in fact linear combinations of a 
function that diverges exponentially and of a function that goes exponentially to 
zero). It is only for a very specific value ofu[{0) that one can make the coefficient 
of the diverging term zero, and have a solution that satisfies mi(oo) = 0. This 
heuristic argument thereby justifies the uniqueness of the number Ui(0) solving 
the problem. 

6.4 Numerical solution 
6.4.1 General principle 

This reformulation of Eqs. (^) and (|37| ) leads to a rather straightforward algo- 
rithm for a numerical solution. Differential problems with two-point boundary 
conditions are usually solved by iterative "shooting" methods where one 
tries to make a guess for the missing initial condition (here, u'i{0)) and then 
corrects this value based on the discrepancy between the resulting end-point 
and the expected one. 

Things are in fact much simpler for an affine equation, since only two trials, 
plus the knowledge of a solution of the complete equation, are enough to find 

-'^The generic solution of 4tu"(t) — ui{t) = is ui (t) = Vi{ciIi{Vt) + C2Ki{Vt)), where 
7i is a modified Bessel function of the first kind. 



16 



the value of u[{0). Indeed, for a differential equation like (|55|), solutions have 
generically the following form: 

ui{t) ^ wit) + aiwi{t) + a2W2it) , (61) 

where w{t) is a particular solution of the full equation, and u'i,2(i) are two 
independent solutions of the corresponding homogeneous equation (i.e. ( |55| ) in 
which one would set uq to zero). Generically, it is convenient to chose these 
functions such that: 

w{0) = , w'{0) = , 
w;i(0) = l, w[{0)^0, 

W2(0) = 0, w^(0) = l. (62) 

If these solutions are known (at least numerically), the condition ui(0) = — 1 
implies 

ai = -1 . (63) 
Then, using wi(+oo) — 0, the second coefficient a2 is given by 

wi(t) — w(t) 

t~, + cx> W2{t) 

and thanks to our choice of initial conditions for w,wi and W2, the value of 
Ui{0) that solves the problem is simply given by 

u[{0) = a2. (65) 

6.4.2 Realistic algorithm 

This algorithm is however not directly applicable in the case of Eq. (^5|), because 
the point t = is a singular point of the equation, and cannot be used to set 
initial conditions. Therefore, we have to chose another point in order to set the 
initial conditions. Let us call this point to > 0, and assume 

wi(to) = l, w[(to) = 0, 

W2{to)^0, w^(to) = l, (66) 



instead of Eqs. (62). From these initial conditions, one must evolve numerically 
the functions w, wi and W2 both forward and backward. Having in mind what 
has been said at the end of section |6.3| , the three functions are going to diverge 
when t +00 (unless one has been very unlucky when choosing the initial 
conditions). The condition ui(+c») = then implies 

= 1 + am + asra , (67) 
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where = Muit^+oa Wi{t)/w(t). This gives a first hnear relation between the 
two cocfBcicnts ai and ck2- Then, from the condition ui(0) = —1, one can 
obtain 

a, = hn. ^Mt)-rAl + wit)) 
t-*o+ riW2(t) - r2Wi[t) 

At this point, the two coefficients ai,2 are known, and it is easy to find 

u'i(O) = lim w'{t) + aiw[{t) + a2w'2{t) . (69) 

This last step does not require any additional calculation, since the derivatives 
of u',wi^2 are known numerically at this point. In summary, this algorithm re- 
duces the problem of solving the integral equation ( |37| ) and then calculating the 
integral over hi Eq. ( p6| ) to the numerical solution of a differential equation 
with three different initial conditions. A numerical analysis of Eq. (|3^) based 
along these lines will be presented elsewhere [p2[. 



7 Conclusions 

In this paper, we have derived a simple sum rule that enables to perform ana- 
lytically some of the integrals involved in the thermal 2-loop photon production 
rate. Several applications of this sum rule have been presented. This sum rule 
also plays a role in the calculation of the collision integral that appears in the 
resummation of the ladder diagrams involved in the calculation of the LPM 
effect. 
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A Asymptotic behavior of F{x) 
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A.l Large x behavior of F{x) 

At large values of the argument x, the asymptotic value of the function F{x) 
introduced in Eq. ( |20| ) is very easy to obtain|^: 

N , tanh~^(u) , 1 

F{x) = e / du ^ — with e = < 1 

Jo W + e X - 1 

^ tanh^^(M)— w u 
du he/ du—::, 







e [1 - ln(2) + 0{e ln(l/e))] + e In (i±^) , (71) 



I.e. 



x~^+oo 2x X \ X J 

Let us add that if the variable x goes to — oo, one can obtain the correct asymp- 
totic behavior by replacing x by —x in the previous expression (in particular 
ln(a;) by In |a;| — ztt) while dropping the imaginary part, so that the previous 
asymptotic formula simply becomes: 

x^-aa 2x X \ X J 

A. 2 Small x behavior of F{x) 

At small values of x, determining the expansion of F{x) requires a little more 
work. Denoting y = l/y/l — x « 1 + a;/2, we have 

1 



F{x) 



2 Jo yi-u 



u/ — u^ 



^ duXiii 



4./n VI 



1 1 



y + u y — u 



(74) 



Let us notice first that the term in l/{y + u) is finite in the limit y 1. 
Therefore, since we do not want to go beyond the constant terms, we can simply 
replace y by 1 in this term and get 

Vl-u/l + u 4/ 1 + v 48 ^ ^ 



^•^By subtracting more of the Taylor expansion of tanh ^(u), one can go one order further 
in this asymptotic expansion, and obtain 

F(x) ~ ^^^^"^ + + ^'^^^'^ + ^~^ +o( ^^^^^ ^ (70) 

x^oo 2x X 3x'^ y J.3 J 
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The other term can be separated in two parts, the first one being 



— I du 



ln(l + u) 

y-u 



y f ln(l + u) - ln(2) ^ y f ln(2) 

4 



y — u 



y-u 



1 f dv , ln(2) , /2\ 



1/2 1/2 

ln(l — f ) /■ ln(l — v) 

4 1 - w J A v{l - v) 



ln(2)^^/2 



^ ln2(2) TT^ ln(2) /2 
up to terms that vanish when x O'^ . Seemingly, the second piece is: 



du 



ln(l - u) 
y-u 



du 



\n{y-u)+\n{l-^) 



y~u 



(76) 



n = l ^ n 

ihl^f^)+^ 

8 \xJ 24 



du 



Collecting all the bits and pieces, we finally obtain: 



F{x) 



1 



In" 



2:^0+ 



)+^ + 0(xln2(l/x)) 



(77) 



(78) 



When a; approaches by negative values, it is sufficient to replace a; by — a; in 
the previous formula (i.e. ln(a;) by In |x| — in) and retain only the real part. 
Since the logarithm is squared, the —iir modifies the constant term as follows: 



F{x) 



I 



■In" 



12 



(xln^(l/a;)) . 



(79) 



B For real photons: Jt = —Jl if M^o = TUg 

In the production of real photons, one can notice numerically that the coeffi- 
cients and are equal if the gluon plasmon mass mg is equal to the quark 
asymptotic mass Moo (for real photons Afcff = Moo). We present here an 
analytic proof of this statement. When Moo — "ig and — 0, we have: 

J,+J, = 7t[F{4/3)-2F(4)] 

L f. tanh-^^ 
= TT 3 / du 2— T 

. 



2 J du 





tanh ^ (u) 



3^2 



(80) 
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Noticing that the function 

G{x) 

obeys the differential equation 



^ , tanh"Vu) 
au- 



(81) 



In 



and solving it, one can prove the following formula 
tanh"^(w) 



(82) 



= ln(2)tan" 



.•tan"^(l/2:) 



r , tanh-Vw) , ,/l\ 1 , Wl+w^) 



d6'ln(cos(6')) . (83) 



Using this intermediate result, it is easy to rewrite + as 



3 / d6lln(cos(6l)) - 2 / rf6'ln(cos(6i)) - 
Jo Jo 

Making then use of the following relations 



(84) 



Tr/2 

de\n{cos{e)) = 


7rln(2) 
2 


7r/3 

d6lln(cos(6')) = 


7rln(2) 
2 


7r/6 

d6'ln(cos(6')) = 


7rln(2) 
2 


7r/6 





7r/3 



d6'ln(sin(6l)) , 



/ dein{sm{e)) + dein{cos{e)) = - rf6'ln(sin(6»)) 



7rln(2) 



it is straightforward to check that 



J+J^ =0 



(85) 
(86) 



when Moo = "^^g- This is an interesting non trivial (and exact) identity which 
cannot be obtained by any simple other method, mainly because there is no 
useful asymptotic formula near the point M^o = TUg. Even if this identity is 
purely anecdotical, its derivation illustrates the power of the sum rule obtained 
in this paper. 
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